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The extremal coefficient function (EOF) of a max-stable process X on some index set T assigns 
to each finite subset A CT the effective number of independent random variables among the 
collection {Xt}t€A. We introduce the class of Tawn-Molchanov processes that is in a 1:1 cor¬ 
respondence with the class of ECFs, thus also proving a complete characterization of the ECF 
in terms of negative dehniteness. The corresponding Tawn-Molchanov process turns out to be 
exceptional among all max-stable processes sharing the same ECF in that its dependency set 
is maximal w.r.t. inclusion. This entails sharp lower bounds for the finite dimensional distribu¬ 
tions of arbitrary max-stable processes in terms of its ECF. A spectral representation of the 
Tawn-Molchanov process and stochastic continuity are discussed. We also show how to build 
new valid ECFs from given ECFs by means of Bernstein functions. 
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1. Introduction 

Besides the class of square integrable processes, the class of temporal or spatial max- 
stable processes is of common interest in stochastics and statistics, cf. [3, 4, 9, 13, 20, 31], 
for example. In spite of considerable differences between these two classes, for example, 
the non-existence of the first moments in case of max-stable processes with unit Frechet 
marginals, connections between the two classes have been made for instance, by the ex¬ 
tremal Gaussian process [25] and the Brown-Resnick process [16] that are parameterized 
by a correlation function and a variogram, respectively. 

Naturally, extremal dependence measures such as the extremal coefficients [26, 28], the 
(upper) tail dependence coefficients [1, 5, 8, 11] or other special cases of the extremogram 
[8] are appropriate summary statistics for max-stable processes. In this article, we capture 
the full set of extremal coefficients of a max-stable process X = {Xt}t^T on some space T 
in the so-called extremal coefficient function (ECF) 9, which assigns to each finite subset 
A of r the effective number of independent variables among the collection {Xt}t^A- We 


This is an electronic reprint of the original article published by the ISI/BS in Bernoulli, 
2015, Vol. 21, No. 1, 276-302. This reprint differs from the original in pagination and 
typographic detail. 


1350-7265 © 2015 ISI/BS 



2 


K. Strokorb and M. Schlather 


introduce a subclass of max-stable processes that is parameterized by the ECF, and thus 
reveal some analogies to Gaussian processes and positive definite functions as follows: 

Among (zero mean) square integrable processes, the subclass of Gaussian processes 
takes a unique role, since it is in a 1-1 correspondence with the set of covariance functions, 
which are precisely the positive definite functions. This fact can be proven by means of 
Kolmogorov’s extension theorem and is illustrated in the following graph: 



In case T is a metric space, the Gaussian process Z* (C) is continuous in the mean square 
sense (and then also stochastically continuous) if and only if the covariance function C 
is continuous if and only if C is continuous on the diagonal (cf. [23], Theorem 5.3.3). 
Well-known operations on the set of positive definite functions C, and hence on the 
corresponding Gaussian processes Z*{C), include convex combinations and pointwise 
limits. Moreover, Bernstein functions play an important role for the construction of 
positive definite functions. 

In our case, the crucial role of zero mean Gaussian processes is taken by the class of 
Tawn-Molchanov processes (TM processes), which are in fact the spatial generalization 
of the multivariate max-linear model of [26]. Using Kolmogorov’s extension theorem, we 
shall see that each ECF 0 (of some max-stable process) uniquely determines a TM process 
X*{9) having the same ECF (Theorem 8). Alongside, we generalize a multivariate result 
[19], Corollary 1, to the spatial setting, proving that the ECFs coincide with the functions 
0 on F{T) (the set of finite subsets of T) that are normalized to 0(0) = 0 and 0({t}) = 1 
for t GT and that are negative definite (or equivalently completely alternating) in a sense 
to be explained below (cf. Definition 6). This can be illustrated in analogy to the above 
sketch: 



Having identified the ECF 0 as a negative definite quantity allows for several immediate 
consequences: First, we obtain an integral representation of 0 as a mixture of maps 
A lyinQ/0 (Corollary 16) and derive a spectral representation for the corresponding 
TM process X*(0) (Theorem 18). Second, we consider operations on ECFs that allow 
to build new ECFs from given ones. We find that ECFs allow for convex combinations 
and pointwise limits (Corollaries 13 and 14) and that the class of Bernstein functions 
operates on ECFs (Corollary 20). We also recover the “triangle inequalities” for 0 from 
[7], Proposition 4, and see that the inequalities therein correspond to three specific choices 
of a Bernstein function, whereas we may plug in arbitrary Bernstein functions to obtain 
further “triangle inequalities” (Corollary 21). 
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Figure 1. Examples of dependency sets in a trivariate setting: a “typical” dependency set K. 
(left) and a dependency set 1C* stemming from a TM process (right). It is shown that 1C C 1C* 
(middle). For further details, see the introduction, Example 31, Lemma 30 and Theorem 32. 


For T being a metric space, we discuss stochastic continuity. The TM process X*{9) 
is stochastically continuous if and only if 0 is continuous (cf. Definition 22) if and only if 
the bivariate map (s,t) e-)- 0({s,t}) is continuous if and only if the bivariate map {s,t) i—>■ 
0({s,t}) is continuous on the diagonal (Theorem 25). 

Finally, we address the exceptional role of the TM processes among simple max-stable 
processes. To this end, Molchanov’s dependency set 1C [19] is transferred to max-stable 
processes X. It comprises the finite dimensional distributions (f.d.d.) of X (Lemma 29). 
Now, let 1C* (9) denote the dependency set of the process X*{9). Then we identify 1C* (9) 
as intersection of halfspaces that are directly given by the ECF 9 (Theorem 32). It turns 
out that 1C* (9) is exceptional among the dependency sets 1C of all max-stable processes 
sharing the same ECF 9, since 1C* {9) is maximal w.r.t. inclusion as illustrated in Figure 1. 
Since inclusion of dependency sets corresponds to stochastic ordering, this observation 
leads to sharp inequalities for the f.d.d. of max-stable processes in terms of its ECF 9 
(Corollary 33). 

The text is structured as follows. After the introductory Section 2, the characteriza¬ 
tion of ECFs and the existence of TM processes is established in Section 3. Section 4 
collects several immediate consequences and related results, while Section 5 exhibits the 
exceptional role of TM processes. Sections 4 and 5 can be read independently. 


2. Foundations and definitions 

2.1. Notation for max-stable processes and ECFs 

A stochastic process X = {Xt}t£T on an arbitrary index set T is said to be max-stable 
if for each n £ N and independent copies X^^\ ..., of X the process of the maxima 
{Vr=i has the same law as {an(t)Xt + &n(t)}teT for suitable norming functions 

a„(<) > 0 and bn(t) on T. Without loss of generality, we shall deal with max-stable 
processes that have standard Frechet marginals, that is, P(At < x) = for t £ T 
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and a; > 0, and set an(t) = n and bn{t) = 0. Such processes are called simple max-stable 
processes. 

It has been shown (cf. [9, 15, 29]) that (simple) max-stable processes X = {Xt}t^T 
allow for a spectral representation (17,^, z/, 14): there exists a (sufficiently rich) measure 
space and measurable functions Vt : —>■ R+ (with f^Vt(uj)i^(duj) = 1 for each 

t G T), such that the law oi X = {Xt}teT equals the law of 

I V UVtioj)] . (1) 

J teT 

Here H denotes a Poisson point process on 1R+ x H with intensity u~‘^du x ^{duj). The 
functions {VtjteT are called spectral functions and the measure v is called spectral mea¬ 
sure. 

In order to describe the finite dimensional distributions (f.d.d.) of X, we shall fix 
some convenient notation first: Let M G T he some non-empty finite subset of T. By 
(resp. [0,oo]^) we denote the space of real-valued (resp. [0, ooj-valued) functions on 
M. Elements of these spaces are denoted by a; = {xt)teM where Xt = x{t). The standard 
scalar product is given through (x,y) = For a subset L C M, we write 1l 

for the indicator function of L in (resp. [0,00]-^), such that {IjtjjtgM forms an 
orthonormal basis of . In this sense, the origin of equals I 0 being zero everywhere 
on M. Using this notation, we emphasize the fact that a multivariate distribution of a 
stochastic process is not any jMj-variate distribution, but it is bound to certain points 
(forming the set M) in the space T. Finally, we consider some reference norm jj • jj on 
(not necessarily the one from the scalar product) and denote the positive unit sphere 
Sm ■= {aG [0,oo)^ : ||a|| = I}. 

In terms of a spectral representation {Q,A,i',V), the f.d.d. of X are given through 

-log¥{Xt<XtA^ M)= I (\J -5^^^V(dw) (2) 

for X G [0, 00 )'^ \ {I 0 }. Alternatively, the f.d.d. of X for a finite subset 0 ^ M gT may 
be described by means of one of the following three quantities that are all equivalent to 
the knowledge of the f.d.d.: 

• the (finite dimensional) spectral measure Hm (cf. [10, 22]), that is, the Radon mea¬ 
sure on Sm such that for x G [0, 00 )^ \ {I 0 } 

-\og¥{Xt<xt,t&M)= [ ( y ^^HMida), (3) 


• the stable tail dependence function Im (cf. [1]), that is, the function on [0,oo)^ 
defined through 


' Sm 


y atXt\HM{da), 


teM 


^m{x) := -logP(Xt < IfxtA^ M) 


( 4 ) 
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• the (finite dimensional) dependency set K-m (cf. [19]), that is, the largest compact 
convex set ICm C [0,oo)^ satisfying 


^M(a;) =sup{(x,?/): 2/e/C m} Vxe[0,oo) 


M 


(5) 


In order to be a valid finite dimensional spectral measure of a simple max-stable random 
vector {Xt}teM, the only constraint that a Radon measure Hm on Sm has to satisfy is 
that 



atHM{da) = 1 


for each t G M. This ensures standard Frechet marginals. 

Given a simple max-stable process X on T, we may assign to a non-empty finite subset 
Act the extremal coefficient 0{A) (cf. [26, 28]), that is 


e{A):= lim 

X—^■oo 


logP(X( < x) 




Hmida) = £m(1a)- 


( 6 ) 


Indeed, the expression logP(Vtg^-^t < x)/ \ogT{Xt < x) does not depend on x and 
equals the right-hand side (r.h.s.) for A C M. Observe that 0{A) takes values in the 
interval [1, j^j], where the value 1 corresponds to full dependence of the random vari¬ 
ables {Xt}t£A and the value jAj (number of elements in A) corresponds to independence. 
Roughly speaking, the extremal coefficient 0{A) detects the effective number of indepen¬ 
dent variables among the random variables {Xt}t^A- It is coherent to set O{0) := 0 to 
obtain a function 0 on T{T), the set of finite subsets ofT. We call the function 


6»: J-(T)^ [0,oo) 

extremal coefficient function (ECF) of the process X. The set of all ECFs of simple 
max-stable processes on a set T will be denoted by 


0(T) = {9:iF{T) -A [0,oo) :0 is an ECF of a simple max-stable process on T}. (7) 


Example 1. The simplest ECFs are the functions 6{A) = JA] corresponding to a process 
of independent random variables, and the indicator function 9{A) = '^A^^es corresponding 
to a process of identical random variables. 


Rather sophisticated examples of ECFs can be obtained using spectral representations 
(O,^, z/, V) of processes X. In these terms the ECF 0 of a process X is given by 


e{A) 


V ^t(w)^i^(dw) 


( 8 ) 


for A G E{T) \ { 0 } and 0(0) = 0. 
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Example 2 (Mixed Moving Maxima - M3 process). Consider the simple max- 
stable stationary process X on that is given through the following spectral represen¬ 
tation 

• {V,,A,i') = (¥xR‘^,iF^B{R‘^),p^dz),wheie d^:) denotes the Lebesgue- 

measure on the Borel-cr-algebra of and where (F, E, p) denotes a measure space 
of non-negative measurable functions on with f{z) dz)p{df) = 1, 

• Vt{{f,z)) = f{t-z) fortsR^ 

then we call X a Mixed Moving Maxima process (MS process) (cf. also [17, 25, 29, 30]). 
Because of (8) the ECF 0 of a Mixed Moving Maxima process X can be computed as 

^{A)= [ [ (y f{t-z))dzp{df) 

for A € E(R‘^) \ { 0 } and 0(0) = 0. In case p is a point mass at an indicator function / 
the bivariate coefficient d({s, t}) will be given by d({s, t}) = 2 — f * f{s — t), where f * f 
means the convolution of / with / and f(t) = f(—t). 

Example 3 (Brown Resnick process). Consider the simple max-stable stationary 
process X on that is given through the following spectral representation (^l.A^v^V): 

• (fl, M, y) denotes the probability space of a Gaussian process W on R'^ with station¬ 
ary increments and variogram 7(2) = E(M4 — Wo)^ for z G R*^. 

• Vt{uj) = exp(Wt(Lo) — a'^(t)l2) for t G R^^, where o'^(t) denotes the variance of W*, 

then we call X a Brown-Resnick process (cf. [16]). Because of (8) the ECF 0 of a Brown- 
Resnick process X is 

9iA)=Eexp(y Wt-a^t)(2^ 

\eA ^ 

for A G J'(R'^) \ { 0 } and 0(0) = 0. Since the f.d.d. of X only depend on the variogram 
7, the extremal coefficient 0(A) will also depend only on the values {7(5 — t)}s,tgA- In 
particular, we have 0({s,t}) = 1 + erf(\/j(s — t)/8) for the bivariate coefficient 0({s,t}), 
where erf(x) = 2)yjn e~* dt denotes the error function (cf. [16]). In case the variogram 
equals 7 (z) = AJjzlP for some A > 0, explicit expressions for multivariate coefficients of 
higher orders up to d -I- 1 can be found in [12]. 

2.2. A consistent max-linear model 

A multivariate simple max-stable distribution is called max-linear (or spectrally discrete) 
if it arises as the distribution of a random vector X of the following form 
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where Z = {Zj}j^-^ is a vector of i.i.d. unit Frechet random variables and where {aij}pxq 
is a matrix of non-negative entries with ~ 1 each row i = 1,... ,p. This is 

equivalent to requiring the spectral measure Hm from (3) for M = {1,....... ,p} to be 

the following discrete measure on Sm 

i=i 

where aj denote the column vectors of the matrix {aij}pxq- Conversely, any discrete 
spectral measure of a simple max-stable random vector gives rise to such a matrix. 
Surely, the ECF of such a random vector X = {Xi}i^M is 

= XI V “d 

j—1 iGA 

for 0 ^ A C M and 0(0) = 0 (cf. (6)). 

In [26], the authors introduce a max-linear model for X* = {X*}i^M where the column 
index j ranges over all non-empty subsets L of M and where non-negative coefficients 
tl are given for each column 0 ^ L G M, more precisely 

= V i e M, with ai,L = TLljgL, 

which is equivalent to 

x:= V tlZl, i&M. (10) 

iGLizM 

The model (10) is simple if and only if = '^lzm ■ izL = f for each i G M. 

It follows from (9) that the ECF of model (10) is 

= X 

LCM : AnL^0 

for 0 ^ Ac XI and O{0) = 0. Now, the interesting aspect of this model (10) with given 
coefficients is that such models are in 1-1 correspondence with ECFs 9 on the finite 
set M (cf. [26], Theorems 3 and 4). Alongside, this leads to a set of inequalities which 
fully characterizes the set of ECFs 0(M) for finite sets M (cf. [26], Corollary 5). An 
alternative proof for these inequalities is offered in [19], Corollary 1, and it is noticed 
therein that they are equivalent to a property called complete alternation (see below). 

As we seek a spatial generalization of these results, let us consider a max-stable pro¬ 
cesses X* = {X^}tzT on an arbitrary index set T, whose f.d.d. for a finite set M are 
precisely of the above form (10), where the coefficients tl now additionally depend on 
M. That means we set the spectral measure Hl^ of the random vector {X^}tzM 

0/LCM 


( 11 ) 
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Figure 2. Spectral measure representation of {X*}t^M for M = {1,2,3} if we choose the ref¬ 
erence norm on to be the maximum norm. In this case, the spectral measure simplifies to 
a sum of weighted point masses on the vertices of a cube: ■ 


such that the f.d.d. of the process X* are given by (cf. (3)) 

-\ognx; <xut€M)= Yl (12) 

0^LCM tGL 

Here M ranges over all non-empty finite subsets of T, which we express as M G X{T) \ 
{ 0 }. Figure 2 illustrates this spectral measure for a trivariate distribution where M = 
{1,2,3} in case the reference norm is the maximum norm. 

Lemma 4. Let T be an arbitrary set and let coefficients rff be given for M G X{T) \ { 0 } 
and L G X{M) \ { 0 }, such that 

(i) > 0 for all M G J-{T) \ { 0 } and L G X{M) \ { 0 }, 

(ii) rff = -|- for all M G XlT) \ { 0 } and L G X{M) \ { 0 } and t G 

T\M, 

(hi) = 1 for all t gT. 

Then the spectral measures from (11) define a simple max-stable pro¬ 

cess X* = {Xf}t^T on T with f.d.d. as in (12). 

Proof. Condition (i) ensures that each spectral measure defines a max-stable dis¬ 
tribution with Frechet marginals. Subsequently, condition (ii) ensures consistency of 
these distributions (i.e., the conditions for Kolmogorov’s extension theorem are satis¬ 
fied). Hence the spectral measures define a max-stable process X* on T. Finally, 
condition (iii) ensures that the process X* has standard Frechet marginals. □ 
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Remark 5. Condition (ii) is equivalent to 

rK= yMeT{T)\{0},0^KcAcM. (13) 

JCM\A 

3. The TM process and negative definiteness of ECFs 

For the following characterization of the set of ECFs 0(F), we use the fact that iF(T), 
the set of finite subsets of T, forms a semigroup with respect to the union operation U 
and with neutral element the empty set 0. The following notation is adopted from [18] 
and [2]. For a function f :T{T) -A R and elements K,L G T(T), we set 

{AKf)iL):=f{L)-f{LUK). 

Definition 6 (negative definiteness/ complete alternation). A functionR(T) -A 
R is called negative definite (in the semigroup sense) on T{T) if for all n > 2, 
{Ki,Kn} C F'(T) and {oi,... ,a„} C R with J2]j=i = 0 

n n 

EE ajOk/iKjU Kk) < 0 . 

1=1 fc=i 

A function ib'.TiT) —>■ R is called completely alternating on FiT) if for all n> 1, 
{Ki, ...,Kr,}C T{T) and K G T{T) 


{AkAk,---Ak^/){K)= Y. {-^P/iKu\jKA<0. 

/C{l,...,n} ^ iel ' 

Because the semigroup (F(r),U,0) is idempotent, these two terms coincide. That 
means ■!/): F'(T) —i R is completely alternating if and only if / is negative definite (in the 
semigroup sense), cf. [2], 4.4.16. 

Example 7 ([18], page 52). An important example of a negative definite (completely 
alternating) function on J-{T) is the capacity functional C :J-{T) —>■ R of a binary process 
Y = {YtjtgT with values in {0,1}, which is given by C(0) = 0 and 

C{A) = P(3t G A such that Yt = 1). 

Now, we are in position to characterize the set 0(T) of possible ECFs on F(T) and 
to define a corresponding max-linear process X*. 

Theorem 8. 

(a) The function 9\E(T) —i R Is the ECF of a simple max-stable process on T if and 
only if the following conditions are satisfied: 
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(i) 0 is negative definite, 

(ii) 0(0) =0, 

(iii) 9{{t}) = 1 for all t €T. 

(b) If these conditions are satisfied, the following choice of coefficients 

rf := \ L) = ^ (-l)l^l+i0((M \L)UI) 

ICL 

VM G T{T) \ { 0 }, 0 ^L = {ti,...,ti}cM 

for model (11) defines a simple max-stable process X* on T which realizes 9 as its 
own ECF 9*. 


Definition 9 (Tawn Molchanov process (TM process) ). Referring to the previous 
work in [6, 19, 26], we will call the simple max-stable process X* from Theorem iS(b) 
Tawn-Molchanov process (TM process) henceforth. 


Proof of Theorem 8. If 0 is an ECF of a simple max-stable process X on T, then 
necessarily 0(0) = 0 and 0({t}) = 1 for alH G T (cf. (6)). Further, it is an application of 
I’Hopitals rule that for A C J-{T) \ {0} 


9{A) = lim 

X—>^oo 


-logP(Vtg^At<a;) 

-logP(At < x) 


lim 


l-P(V,g^A«<cr) 

1 - P(At < x) 


P(3t G A such that Xt > x) 

hm - — -;- 

c-s-oo P(^t > x) 


lim 




(14) 


where denotes the capacity functional for the binary process Yt = lxt>x and p^^'> = 
WYt = 1 — Since negative definiteness respects scaling and pointwise limits, negative 

definiteness of 0 follows from Example 7. This shows the necessity of (i)-(iii). 

Now, let 0 :iF{T) —>■ M be a function satisfying conditions (i)-(iii) and let the coefficients 
be given as above. We need to check that they fulfill the (in)equalities from Lemma 
4. Indeed we have: 

• The inequalities = — Ajj^} • • • A|(j}0(M\L) > 0 follow directly from the complete 

alternation of 0 that is equivalent to (i). 

• From the definition of Aj^j we observe 


= -A|*} A|*,}... A|*,}0((M U {t}) \ {L U {t})) 

= -A|t,} • • • A|t,}0(M \ L) + Ap,} ■ • • App0(M U {t} \ L) 

_ .^M MU{t} 

— ^L-^L 


• For t G T, we have = 0({t}) = 1 because of (iii). 
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Thus, the coefficients define a simple max-stable process X* on T as given by model 
(11). Finally, we compute the ECF 9* of X* and see that it coincides with 9: For the empty 
set, we have 9*{0) = 0 = 9{0) because of (ii); otherwise we compute for A C iF{T) \ {0} 
that 

( 6 )^ 11 ) Y. Tt= Y. 

0^LcA 0^LCAICL 

= E W E E 9iK)i-{-lK=A))=0iA). 

0^KCA 0^LCA 0^KCA 

A\LCK 

This shows sufficiency of (i)-(iii) and part (b). □ 

Theorem 8 is in analogy to the following standard result for Gaussian processes (as 
illustrated in the sketches in the Introduction): 

(a) A function C :T x T ^ M. is a covariance function if and only if it is positive 
definite. 

(b) If C: T X T —>■ R is positive definite, we may choose a (zero mean) Gaussian process 
which realizes C as its own covariance function. 

Both statements are intrinsically tied together. When proving them by means of Kol¬ 
mogorov’s extension theorem, one proceeds in the same manner as we did for Theorem 
8. The necessity of positive definiteness of covariance functions is easily derived even for 
the bigger class of square-integrable processes, whilst sufficiency can be established by 
showing that Gaussian processes can realize any positive definite function as covariance 
function. In some points (such as continuity relations), this analogy will be deepened. 
Other aspects (such as the exceptional role of dependency sets in Section 5) seem un¬ 
suitable for a direct comparison. 

Remark 10. In order to incorporate stationarity w.r.t. some group G acting on T (for 
example, R'’* acting on R'^ by translation), we just have to add the following condition 
(iv) 9{gA) = 9{A) for all A € .F(T) \ {0} and for all g GG. Then the process X* will be 
stationary w.r.t. this group action. 

Remark 11. Instead of requiring the max-stable processes in Theorem 8 to have stan¬ 
dard Frechet marginals everywhere, we can admit a different scale at different locations, 
that is, F{Xt < a;) = exp(—st/a;) for a positive scaling parameter st for t G T. In that case 
Theorem 8 holds true without condition (iii) and the word “simple”. To make sense of the 
EGF as in (6) in this case, either use a reference point tGT oi set logP(At < x) = —1/x 
in the denominator. Beware of that the EGF 9 cannot be interpreted as the number of 
independent variables anymore in this case. 

Remark 12. In [26] , the last issue of the proof is derived for finite sets T by a Moebius 
inversion. The relation to the proof therein becomes more transparent if we compute 
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9* (A) for A C M from the coefficients for arbitrary M Z) A instead of 

M = A: 

E E E ^Kuj= E ^L- (15) 

0^KcA 0^KcAj(ZM\A LCM: LnA^0 


4. Direct consequences of Theorem 8 

Here, we collect some direct consequences of the above Theorem 8. Therefore, note that 
the first part of Theorem 8 can also be expressed as (cf. (7)) 

0(T) = {0:T(T) -A [0, oo): 6 is negative definite, 0(0) = O,0({f}) = 1 for t € T}. (16) 


4.1. Convexity and compactness 

Corollary 13. The set of ECFs 0(T) is convex. 

Proof. This can be seen directly from (16) since all involved properties are compatible 
with convex combinations. As a constructive argument, use the fact that the ECF of the 
max-combination aX V (1 — a)Y of two independent simple max-stable processes X and 
Y on T is the convex combination of their ECFs for a G (0,1). □ 

Corollary 14. The set of ECFs 0(T) is compact w.r.t. the topology of pointwise con¬ 
vergence. 

Proof. The topology of pointwise convergence on is the product topology. Since 

0(0) = 0 and 0(A) G [1, |A|] for 0 G 0(T) and A G E(T) \ {0}, the set 0(T) is a subset 
of the product space 

{o}x n 

AeJ^(r)\{0} 

which is compact by Tychonoff’s theorem. Moreover, since elements of 0(T) are com¬ 
pletely characterized by finite dimensional equalities and inequalities involving < only 
(stemming from (16)), the set 0(T) is closed. Hence, 0(T) is compact. □ 

Remark 15. Note that even though we say “the topology of pointwise convergence”, the 
“points” meant here are indeed elements of J'(T), that is, finite subsets of T. In particular 
it follows from the compactness of 0(T) that 0(T) is sequentially closed. That means if 
(0ri)nGN IS & Sequence of ECFs such that 0„(A) converges to some value f{A) for each 
A G R{T), then / is an ECF. 
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Another consequence of Theorem 8 is that ECFs allow for an integral representation as 
a mixture of functions A i—>■ l^nQ/ 0 ) where Q is from the power set of T. To be more 
precise, let us denote the power set of T by V{T) and consider the topology on V{T) 
that is generated by the maps Q ^ lAnQ /0 for A € iF(T) or equivalently (since iF(T) 
is generated by the singletons {{t}}tgT) the topology on V{T) that is generated by the 
maps Q ItgQ for t GT. Identifying 'P{T) with {0,1}^, this space is also known as 
Cantor cube. As in [2], Definition 2.1.1, a measure p, on the Borel-cr-algebra of V{T) 
w.r.t. this topology will be called Radon measure if p is finite on compact sets and p is 
inner regular. 

Corollary 16. Let 6 G Q(T) be an ECF. Then 6 uniquely determines a positive Radon 
measure p on V(T) \ {0} such that 

0{A)=p({QGViT)\{0}-.AnQ^0})= [ lAnQ^0lJ-idQ), 

where = 1 for t GT . The function 6 is bounded if and only if p{V{T) \ {0}) < oo. 

Proof. Since 0 is negative definite (Theorem 8 ) and J-{T) is idempotent, we may apply 
[2], Proposition 4.4.17. It says that 0 uniquely determines a positive Radon measure p on 
J-{T) \ {1}, where iF{T) denotes the dual semigroup of iF{T) (cf. [2], 4.2.1 and 4.4.16), 
such that 0{A) = 'p{{p G .F(T) \ {1} | p{A) = 0}). The function 0 is bounded if and only 
if/I(.^)\{l})<oo. 

Now, it can be easily seen that semicharacters on E{T) are in a 1-1 correspondence 
with subsets of T via J^{T) 9 p —>■ {t G T: p{{t}) = 0} G V{T) and P(T) 9 Q —>■ l(.)nQ =0 G 
J-[T). Here the constant function 1 corresponds to the empty set. Moreover, the topology 
considered on F{T) is the topology of pointwise convergence. Transported to V{T) this 
is the topology generated by the maps Q i -9 lyinQ #0 for 4 G F{T). Let p denote the 
Radon measure p transported to P(T) \ {0}. Then the corollary follows. □ 

Remark 17. In case T = M is finite, we have that V{M) =F{M) carries the discrete 
topology and 

0{A)=p{{QgF{M)\{ 0}:A0Q^0})= T{{Q})^AnQ^0. 

Qe.F(M)\{0} 

A comparison with (15) reveals that p{{Q}) = Tq . In this sense, the coefficients Tq of 
the max-linear model (11) can be interpreted as finite dimensional “Fourier coefficients” 
of the negative definite function 0. 

The integral representation of the ECF 0 also yields a spectral representation for the 
corresponding TM process X*. 
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Theorem 18. The TM process X* = {X^}t£T with ECF 9 has the following spectral 
representation (Jl, At v^V) (cf. (1)): 

• is the measure space lfP{T),B{V{T)), p) from Corollary 16, 

• ^t(Q) = IteQ- 

Proof. We need to check that the f.d.d. of X* satisfy (2). The f.d.d. of X* are given by 
(cf. (12)) 


-\ognx;< Xt,t&M)= Y. 

0^LCM teL 

where the coefficients can be computed from the ECF 6 as in Theorem 8(b) and 9 
satisfies the integral representation from Corollary 16, that is, 

ICL 

= f l((M\L)u/)nQ7^0M(dQ). 

ICL JnT)\{0} 


Using the identity 

E(_l)l^l + ^^((M\£)U7)nO^0 

ICL 

= + l/nQ/0 - l(M\L)nQ#0l/nQ/0) 

ICL 

= 0 • l(M\L)nQ/0 + (1 - l(M\L)nQ/0) E (-l)dl+4/nQ/0 

ICL 

= l(M\L)nQ= 0 lLcQ = lL=MnQ, 


we obtain that 


rf=[ 

JnT)\{0} 

It follows that the f.d.d. of X* satisfy 


lL=MnQM(dQ). 


-log¥{Xf <xt,t€ M) = f Y ^L=MnQ\/—liidQ) 

Jv{T)\{0} 0 ^lcM tcL 



as desired. This finishes the proof. 


□ 
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4.3. Triangle inequalities and operation of Bernstein functions 

In [7], Proposition 4, it is shown that an ECF 9 on iF{T) satisfies the following bivariate 
inequalities for r,s^t€T: 

9{{s,t}) < e{{s,r})9{{r,t}), 

0({s,n)“ < 0{{s,r}r + 0{{r,t}r - 1, 0 < a < 1, 

Oi{s,t}r > 0({s,r}r + e{{r,t}r - 1, a < 0. 

These inequalities have in common, that they are in fact triangle inequalities of the form 

g o r]{{s,t}) < g o r]{{s,r}) g o r]{{rA}), (17) 

if we rewrite them in terms of 77 := 0 — 1 and 

g{x) =log(l + a;), 

g{x) = (1 + xY — 1, 0 < a < 1, 

g{x) = 1 — (1 + xY, a < 0. 

These functions g have in common that they are in fact Bernstein functions. 

Definition 19 (Bernstein function). A function g: [0,oo) —>■ [0,oo) is called a Bern¬ 
stein function if one of the following equivalent conditions is satisfied (cf. [2], 4.4-3 and 
page I 4 I) 

(i) The function g is of the form 

nOO 

g{r)=c-\-br-\- / (1 - e“'^’’);/(dA), 

Jo 

where c,b>0 and n is a positive Radon measure on (0,oo) with Y4_j^(dA) < 
00 . 

(ii) The function g is continuous and g G C''^((0, 00 )) with g >0 and (—> 0 
for all n>0. (Here, denotes the nth derivative of g.) 

(hi) The function g is continuous, g >0 and g is negative definite as a function on 
the semigroup ([0,oo),-|-,0). 

For a comprehensive treatise on Bernstein functions including a table of examples, see 
[24] . Bernstein functions play already an important role in the construction of advanced 
Gaussian processes by generating novel covariance functions from given ones, cf. [32] and 
[21]. Here, we see that they are equally useful for generating new ECFs from given ECFs 
and correspondingly new Tawn-Molchanov processes from given ones. 
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Corollary 20. Let T be a set and 9 £ Q(T) an ECF. Let g be a Bernstein function 
which is not constant. Then the function on T{T) 

g{e{A))-g{f}) 


is again an ECF in Q{T). 

Proof. The result is immediate from Theorem 8, since Bernstein functions operate on 
negative definite kernels (cf. [2], 3.2.9 and 4.4.3). □ 

For instance, if 9 is an ECF, then also log(l + 0)/log(2) or ((0 + a)'^ — a'^)/((l + a)'^ —a^) 
are ECFs for 0 < g < 1 and a > 0. Finally, we show that (17) holds true for arbitrary 
Bernstein functions. In fact, the result of [7], Proposition 4, can be generalized to the 
following extent as a corollary to Theorem 8. 

Corollary 21. Let 9 £ 0(T) be an ECF. Set r] := 9 — 1 and let g be a Bernstein function. 
Then we have for A,B,C£ J-{T) \ {0} that 

g o ri(A U B) < go rjiJC) + go r]{A LI B) < go g^A LI C) + g o g(C U B). 

Proof. Since 9 is an ECF, it is negative definite (cf. Theorem 8). Subtracting 1 does not 
change this property. Notice further that 9 takes values in {0} U [1, oo), where the valne 
0 is only attained for the empty set 0 (the neutral element of F{T)). Thus, the function 
g := 9 — 1 :J^(T)\{0}—J-M is negative definite and takes values only in [0, oo). Applying 
a Bernstein function g does not change this property (cf. [2], 3.2.9 and 4.4.3). By [2], 
8.2.7, this also means that f '■= gog : LF{T) \ {0} —>■ K. is negative definite on LF{T) \ {0}. 
Since we have also / > 0 on J-{T) \ {0}, we may derive for A,B,C £ J-{T) \ {0} 

/(C) + /(A UB)- /(A U C) - f{C U B) 

= {f{C) - /(AU C) - /(CU B) + /(AU i?U C)) + (/(A U B) - f{AU BUC)) 

= AaAb/(C) + Ac/(A U B) < 0 

as desired. This finishes the proof. □ 

4.4. Stochastic continuity 

In this section, we require T to be a metric space. We need to define the notion of 
continuity that we will use in connection with ECFs 9:LF{T) —>■ [0,oo). Therefore, let 
/ :J-{T) —>■ R be a function on the finite snbsets of T. Then / induces a family of functions 
{/*'™^}m>o where :T™ —>■ R is given by 

f^'^\ti,...,tm)= f{{ti,...,tm})- 
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Definition 22. Let f:iF(T) —^ R 6e a function on the finite subsets of a metric space 
T. We say that f is continuous if all induced functions > M are continuous for 

all m > 0, where T™ is endowed with the product topology. 

Lemma 23. Let X = {Xt\t^T be et simple max-stable process with ECF 6. Then the 
following implication holds: 

X is stochastically continuous => 0 is continuous. 

Proof. Stochastic continuity of X means that for any e > 0, for any t gT and sequence 
i(") have P(|X((„) — Xt\ > e) —?► 0. From this, we can easily derive that for any 

e > 0, any m € N, any (ti G T™ and a sequence (4"^, ■ • ■, 4"^) —(<i,..., tm), also 

P(||(X^(n) — XtJlfLiW > e) —>■ 0 for any reference norm || • || on R™. The latter implies the 
corresponding convergence in distribution: F^^(n.) ^(n)^ Since log: 

[ 0 , oo)™ —R is monotone and homogeneous, we have that for x > 0 the point (x,..., x) G 
(0,oo)™ is a continuity point of F(^tx,...,t^) (cf- [22], page 277). Thus, the induced function 
0 ("i) on T™ is continuous, since ... ,tm) = —a:logF'(tj_...^t^)(x,... ,x). Hence, 9 is 

continuous. □ 

Second, we prove the following upper bound that shows that stochastic continuity of 
the TM process X* is indeed controlled by the bivariate extremal coefficients. 

Lemma 24. Let X* = {Xf^t^T be the TM process with ECF 9. Set rj := 9 Then we 
have for any e > 0 

P(|X: - X;| > £) < 2(1 - < ^r,({s,t}). 

Proof. Let e > 0. We will prove the statement for 2e instead of e. Therefore, consider the 
following disjoint events on a corresponding probability space (H, .A, P) for fc = 0 , 1 , 2 ,... 

Hfe := {w G H: (X:(uj), X;(u:)) G (ke, (k + 2)e]^ \ {{k + l)e, (k + 2)ef}. 

The disjoint union IJ^q Ak is a subset of {w G H: |2f*(a;) — Xf{ijj) \ < 2e} and so 

C OO \ 00 n 

IjHfe =^P(Afe) = J^^P(Afc). 
k=0 / k=0 k=0 

From (12) and Theorem 8 , we see that the bivariate distribution of the process X* is 
given by 

-logP(X:<x,X;<y) = ^^ii^ + ^. (18) 

X V 2 / X Ay 

For further calculations, we abbreviate for p, g G N U {0} 

B{p,q)-.= nx:<p-e,x:<q-e). 
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Note that B{p,q) = B{q,p) and B{p,0) = 0. With this notation, we rearrange 

n n 

^ ^ = — B(^ti 1,77- + 1) + 2 ^ + 2,A: + 1) — 5(A: + 2, A:)]. 

k^O k^O 

For the second summand, we have (cf. (18)) 


^ ^ [B{k + 2,A: + 1) — B(^k + 2, k)] 


( 18 ) 


="E 

fc=0 


exp - 


Vi{s,t}) , 1 


/c + 2 /c + 1 


-exp - 






k + 2 



exp - 


k 2 k 


1 


ik + l)ej ke 


>5;^exp(- 




2 e 


exp 


1 


(fc + l)e 


“PI -5 


= expl-%ii) |exp 


2e J (n + l)e J ' 

Finally, 

p(|x;-x;i>2£) 

n 

= 1 - P(|X* - X;| < 2£) < 1 - lim VP(x4fe) 

n—^oo • ^ 

fe =0 

n 

= 1+ lim i3(n + 1, n + 1) — 2 lim 'y^[B{k + 2,k + 1) — B{k + 2,k)] 

n—^oo n—^oo ^ 

fc =0 

< 1 + li,„ expf-5<Wi^) - 2 lim (expf-5<hlh) exp(-^^ 

ra->oo (n+l)£ J n-s-oo^ 2s J \ (n + l)e 

= 2-2exp(-^^t+dl) < ^v{{s,t}). 


This finishes the proof. 


□ 


Theorem 25. Let X* = {X(*}tGT be the TM process with ECF 6 . Then the following 
statements are equivalent: 

(i) X* is stochastically continuous. 

(ii) 6 is continuous. 

(iii) The bivariate map (s, t) i—> 0({s, t}) is continuous. 

(iv) The bivariate map {s,t) i—0({s,t}) is continuous on the diagonal. 
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Proof. The implication (i) (ii) follows from Lemma 23. Clearly, continuity of 9 implies 
continuity of the induced function d^'^^{s,t) ■= which implies continuity of 9^‘^'> 

on the diagonal. This shows the implications (ii) (iii) and (hi) (iv). Finally, the 
implication (iv) (i) follows from Lemma 24, since ri{{t,t}) = ^({0) — 1 = 0. □ 

5. Dependency sets — the special role of TM processes 

In this section, we show that the TM process X* with ECF 6 is exceptional among all 
max-stable processes sharing the same ECF 9 as X* in the sense that its dependency set 
/C* (to be introduced below) is maximal w.r.t. inclusion. 

Therefore, recall that for a finite non-empty subset M CT the dependency set /Cm of 
{Xt}tGM is the largest compact convex set Km C [0,oo)*^ satisfying (cf. (5)) 

^M{x)=snp{{x,y)■.y€KM} VxG[0,oo)^. 

The closed convex set Km may also be described as the following intersection of half 
spaces (cf. [27], Section 1.7): 


Km= f] {y&[0,oo)^ :{x,y)<iM{x)}. (19) 

x£Sm 


Example 26 ([19], Example 1 and Proposition 2). The simplest examples for de¬ 
pendency sets Km are the unit cube [ 0 , 1 ]^ corresponding to a collection of independent 
random variables {Xt}t^M and the cross-polytope := {x G [ 0 ,oo)^ -J^teM^t < 1 } 
corresponding to identical random variables {Xt}teM- Any dependency set Km satisfies 

C Km C [0,1]^. 


Example 27 (Brown Resnick process/Hiislei—Reiss distribution) . The f.d.d. of 
a Brown-Resnick process (cf. Example 3) are the multivariate Hiisler-Reiss distributions 
(cf. [14]). In the bivariate case, when M = {1,2} consists of two points only, the distri¬ 
bution function of a Hiisler-Reiss distributed random vector {Xi,X 2 ), standardized to 
unit Frechet marginals, is 


- 10gP.y(Xl <Xl,X2< X 2 ) 


Xi V 2 1/7 J X 2 \ 2 


log(a:i/a:2) \ 

J 


for xi,a :2 > 0. Here $ denotes the distribution function of the standard normal distri¬ 
bution and the parameter 7 is the value of the variogram between the two points (cf. 
Example 3). Figure 3 illustrates, how the corresponding dependency sets range between 
full dependence (7 = 0 ) and independence (7 = 00 ). 


In order to define a single dependency set for a simple max-stable process comprising 
all multivariate dependency sets, we write 

pr^:[ 0 ,oo)^-;> [ 0 ,oo)“, 


(xt)t^T {xt)teM 
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Figure 3. Nested dependency sets of the bivariate Brown-Resnick (resp. Hiisler-Reiss) 
distribntion where M = {1,2} (cf. Example 27). The dependency sets grow as the parameter 7 
increases. They range between full dependence (7 = 0) and independence (7 = 00 ). 



for the natural projection. 

Definition 28. Let X he a simple max-stable process X = {Xt}t^T o,nd denote for finite 
M G iF{T) \ { 0 } the multivariate dependency set of the random vectors {Xt}t(£M by K-m ■ 
Then we define the dependency set /C C [0,oo)^ of X as 

^■= n 

MgJ^{T)\{0} 

Analogously to (5), the dependency set K. may be characterized as follows. 

Lemma 29. The dependency set K, of a simple max-stable process X = {Xt}teT is the 
largest compact convex set K. C [0,oo)^ satisfying 

^m(x) = supj ^ xtVt : y e /c| Vx € [0, oo)^V0 gT{T), (20) 

where £m is the stable tail dependence function of {Xt}t^M■ 

Proof. Convexity of K. follows from the convexity of each /Cm and from the linearity of 
the projections prM for M £ X{T) \ {0}. Since /C{t} = [0,1] is the unit interval for each 
t £ T, the set AC is contained in the compact space [0,1]^. Moreover, AC is closed as the 
intersection of closed sets, hence AC is compact. 

Next, we prove that ACm = PI'm(^)- By definition of AC it is clear that prM(AC) C ACm 
for M £ iF{T) \ {0}. To prove the reverse inclusion, let j/m be an element of ACm and set 
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ViVM) := pi M {{vm}) r\JC = pr {{pm}) r\ICr\[0, 1]^. We need to show that V{yM) ^ 0- 
Denoting V{yM,A) :=pil^{{yM})r\pT^^{JCA) n [0,1]^, we see that 

V{yM)= n V{yM,A). 

AeT{T)\{0} 

Note that each V{yM,A) is a closed subset of the compact Hausdorff space [0, l]"^. 
Therefore, it suffices to verify the finite intersection property for the system of sets 
{V(vm, A)} a^j^(t)\{ 0 } in order to show that VijjM) ^ 0- But this follows from the con¬ 
sistency of the finite dimensional dependency sets {K-a}a^f(t)\{ 0 } as follows: As [19], 
Section 7, Proposition 8, essentially says, we have that if A and B are non-empty finite 
subsets of T with A<Z B, then 1C a is the projection of /Cs onto the respective coordinate 
space. In particular, prj^^(/CB) C pi'^{1Ca) and prj^^dy^}) fl prj^^(/C_B) fl [0,1]'^ ^ 0 for 
Pa S ICa- Now, let Ai,..., be non-empty finite subsets of T. Then 

0 7^ P^M ({j/m}) n Pr^uudi A,) 1]’^ 

k k 

c p^M ({yM}) n Pi pi^]{lCAi) n [o. 11^ = 0^ {yM,A,), 

2—1 2=1 

as desired and we have shown that /Cm C pIJ^^{IC). Both inclusions give K-m = pinti^)- 
By definition, we have iMix) = sup{(a:, y): y € /Cm} for x € [0, oo)^. Thus, (20) follows 
from /Cm =prM(/C). 

Finally, let £ C [0,oo)^ be also convex compact and satisfying (20) with AC replaced 
by £. Then it follows immediately that pr^(£) = /Cm for any non-empty hnite subset 
M CT. We conclude that £ C AC by definition of AC. This finishes the proof. □ 

In particular, the ECF 0 of a simple max-stable process X = {Xt}teT can be expressed 
in terms of the dependency set AC of A as 

0(A) = sup| ^ It :a: G AC |. (21) 

In order to make statements about the dependency sets AC of processes X = {AtjteT 
in terms of the ECF 0, we introduce the following notation: For any non-empty finite 
subsets A of T, we set the halfspace 

7/^(0) := ji G [0,oo)^ :^a:t < 0(A) 

^ t^A 

that is bounded by the hyperplane 
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Lemma 30. Let K, be the dependency set of a simple max-stable process X = 
with ECF 6. Then the following inclusion holds 

ICC f| UAiO). 

AGCiT)\{0} 

On the other hand for each A & !F{T)\{0} there is at least one point x'^ in the inter¬ 
section 

x^GicnSAie). 

Proof. Let A G X{T) \ {0} and x G 1C. Then the assumption ^ contra¬ 

dicts 9{A) = sup{^(g^xt :xGX}> 9{A) (cf. (21)). So J2t&A — ^(^)- This proves the 
inclusion. Second, since /C is compact and the map [0,oo)^ 9 a; —>■ i® continuous, 

we know that it attains its supremum at some x-^ G X. □ 

Example 31. We give a simple multivariate example for Lemma 30 (as illustrated in 
Figure 1 in the introduction for the trivariate case): The Euclidean norm iMix) = ||a :||2 
is a stable tail dependence function on [ 0 , 00 )-^ (cf. [19], Example 2) and defines a simple 
max-stable distribution (cf. (4)) with ECF 9{A) = for Ac M, such that 

nAi9) = {xG [0, 00 )^ : (x, 1 a) < V\A\}, 

£a(9) = {a: e [0, 00 )“ : (a:, 1 a) = v^} 

for 0 ^ Ac M. It can be easily seen that for x G [0,oo)'^ \ {I 0 } 

iM{x) = ||a :||2 = {x,xl\\x\\ 2 ) = sup{(a;, y): y £ B+}, 

where 5+ := {y G [0,oo)^ : ||y ||2 < 1} denotes the positive part of the (Euclidean) unit 
ball. So, the dependency set K. is clearly in this case. Now, the planes £a(9) are 
tangent to the boundary of with common points x^ = lA/^/\A\ for 0 ^ A C M, 
which makes it easy to see that Lemma 30 holds true in this example. Figure 1 shows the 
dependency set K, = (left) and the intersection of halfspaces bounded by the planes 
£a(,9) (right). In the middle it is illustrated that this intersection contains and the 
points x'^ are marked. 

The following theorem shows that the inclusion from Lemma 30 is sharp and attained 
by TM processes. Figure 4 illustrates the dependency set of a trivariate distribution of a 
TM process. 

Theorem 32. Let X* be the dependency set of the TM process X* = {Xf }teT with ECF 
9. Then 

x*= f| nA{9). 

AeC{T)\{0} 
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Figure 4. Dependency set JC* of the random vector {Xf }t£M for M = {1,2,3}. The dependency 
set /C* is bounded by the hyperplanes £a{Q) that are given by the equations ~ 

where 6 denotes the ECF of X*. The coefficients for L G T{M) \ {0} and t£L turn up as 
lengths of the resulting polytope 1C* (cf. Theorem 8 (b) and Theorem 32). 


Proof. First, we prove the theorem in the case, when T = M is finite and 1C* = IC^: 
Therefore, write 

JCm ■= Pi HAid) = {x € [0,oo)^^ : {x,1a) <0{A) for a\\ 0 ^ Ac M}. 

0^AcM 

The inclusion C £m is proven in Lemma 30. So, it remains to show the other inclusion 
Cm C /C^. Due to (19), we have that 

^M= n {ye [0,oo)^X(x,y) <£^( 0 :)}, 
x£Sm 


where 

4 /(®) = '^L\/xt 

0^LCM teL 

is the stable tail dependence function of here expressed in terms of the coeffi¬ 

cients from Theorem 8 (b) (cf. (12)). Thus, it suffices to show the following implication 
in order to prove Cm C : 

xCSm and y C Cm (x, y) < ^^(a;). 

We now prove this implication: Without loss of generality, we may label the elements 
of M = {ti,..., tm} such that > Xtj > • • • > Xj^. Then we may write x = (xt)t^M S 
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Sm C [ 0 , 00 )''^ as 

X = xt„ 1 m + - a;t„) +-^ {xt2 - 2 ^* 3 ) l{ii.t2} + ( 2^*1 - xt^) 

>0 >0 >0 >0 

Taking the scalar product with y G Cm , we conclude 

{x,y) < xt^6{M) + -xtJ)0{M\{tm}) 

+-H (xt2 - Xt^)e{{txM}) + (a;ti - xt^)e{{ti)) (22) 

= Xt^{9{M) - e{M\{tm})) H- \- Xt^{0{{ti,t2}) - 0i{ti})) + Xt^0{{ti}). 

On the other hand the stable tail dependence function i*j^ is by this ordering of the 
components of x given as 

m 

(*m{.x)= x^\/xt=Y^U 

0^LCM KL i=l 

From (15), we see that this expression coincides with the r.h.s. of (22). Thus, we have 
our desired inequality {x^y) < i^{x). This finishes the proof in the case, when T = M is 
Hnite. 

Otherwise, the definition of the dependency set 1C* and the result for finite M give 

^*= n prM(^M)= n n 

MeJ^(T)\{0} MeT(T)\{0}0=^A(lM 

where 'H^{0)={x€ [0, 00 )^ : ^ ^(^)}- Since {0)) = % a{0) for <Z 

Ag the claim follows. □ 

So, if we fix the EOF 0 of a simple max-stable process on T, then the TM process 
yields a maximal dependency set /C* w.r.t. inclusion, that is 

1C*= U /C. (23) 

K dependency set 
with the same ECF as K* 

Now, inclusion of dependency sets corresponds to stochastic ordering in the following 
sense (cf. [19], page 242): If 1C' and 1C" denote the dependency sets of the simple max- 
stable processes X' and X" respectively, then /C' C /C" implies 


LCM : ti 




rf 


P(X' < xt,t C M) >P(X'' < xt,t C XI) VxG [0,oo)^ 
for all XI G J^(T) \ { 0 }. This leads to the following sharp inequality. 
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Corollary 33. Let X = {Xt}t£T be a simple max-stable process with ECF 6. Let M be 
a non-empty finite subset of T. Then 

P(Xt<xt,teM)>exp^- ^ V Vx e [0,oo)^, (24) 

where the coefficients depend only on 9 and can be computed as in Theorem iS(b). 
Equality holds for the TM process X*. 


Example 34- Let us abbreviate pA '■= d{A) — 1. In the bivariate case, the inequality 
(24) reads as 


P(Xs < Xs,Xt < xt) > exp 
= exp 


Pst 

Xs V Xt 


Xs A Xt 


??st + 1 \ 

Xs A Xt J 


exp 



1 

Xs 



Indeed this inequality is much better then the trivial inequality P(Xs < Xs, Xt < Xt) > 
P(^s < Xs A Xt,Xt < Xs A Xt)-, which can be written in the above terms as 


P(Xs <Xs,Xt< Xt) > exp 


Pst P 1 ^ 
Xs A Xt J 


Further note that pst = ^*({S)i}) ~ 1 can be interpreted as a normalized madogram: 


(14) 

Pst = hm 
x—foo 


P(-^S >XOTXt>x) 
nXt > x) 


-1= lim 




2E1 


Xt>x 


If we additionally take into account that (cf. [26], inequality (13)) 


Prs V Pst V Prt V {prs 4“ Vst 4“ Prt 1) ^ Prst C: {Prs 4“ Pst) A {rjst 4“ Prt) A {j]rt 4“ Prs): 


we obtain from (24) the following (sharp) inequality for the trivariate distribution of a 
simple max-stable random vector (Xr, Xs, Xt) from bivariate quantities: 


P(Xr <Xr,Xs<Xs,Xt<Xt) 

1-prs'^ Pst V Prt 


> exp — 


Pr 


Xr A Xs A Xt 

, 9st 


4" (^Orst A 1) 


1 


1 


1 


Xr AXs Xs A Xt Xr A Xt 


^rt 


^rst 


1 


Xr/\Xs XgAXt XrAXtJ \Xr 

where Orst ■= (Prs EVst) A {prs + Prt) A {pst +Prt)- 


1 

Xs 


1 

Xt 


Pst 

Xr 


Prt 

Xs 


Prs 

Xt 


Thus, if one can handle the ECF of a max-stable process, sharp lower bounds for its 
f.d.d. are available. However, beware that higher variate cases of these inequalities will 
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be numerically unstable. 

Remark 35. It is an open problem and it would be interesting to know whether there 
exist also minimal dependency sets in the sense of (23) and if they would help to better 
understand the classification of all dependency structures. In view of Lemma 30 and 
Theorem 32 a very naive idea would be to take one point from each of the sets /C* fl Ea 
where A £ R{T) \ {0} and then to take the convex hull with 0 included. However, this 
fails to be a dependency set in dimensions |T| > 3, since it is not even a zonoid, which 
would be necessary (cf. [19]). 
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